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Impact of doping on the carrier dynamics in graphene 
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We present a microscopic study on the impact of doping on the carrier dynamics in graphene, in particular 
focusing on its influence on the technologically relevant carrier multiplication in realistic, doped graphene sam¬ 
ples. Treating the time- and momentum-resolved carrier-light, carrier-carrier, and carrier-phonon interactions 
on the same microscopic footing, the appearance of Auger-induced carrier multiplication up to a Fermi level of 
300 meV is revealed. Furthermore, we show that doping favors the so-called hot carrier multiplication occurring 
within one band. Our results are directly compared to recent time-resolved ARPES measurements and exhibit 
an excellent agreement on the temporal evolution of the hot carrier multiplication for n- and p-doped graphene. 
The gained insights shed light on the ultrafast carrier dynamics in realistic, doped graphene samples. 


I. INTRODUCTION 


II. THEORETICAL APPROACH 


A number of theoretical and experimental studies has been 
performed aiming at a thorough understanding of the carrier 
relaxation dynamics in optically excited graphene.*"'^ Most 
of these studies focus on the ultrafast Coulomb- and phonon- 
induced carrier dynamics without considering the influence 
of doping in the investigated graphene samples. A non-zero 
Fermi level can have a crucial impact on the relaxation dy¬ 
namics via a significant increase of the scattering phase space 
and via the enhancement of Pauli blocking. A first experi¬ 
mental time-resolved ARPES study has been performed ad¬ 
dressing the doping dependence of carrier multiplication of 
graphene.*^ The underlying elementary processes determin¬ 
ing the observed different behavior for p- and n-doped sam¬ 
ples have not been microscopically investigated, yet. 

In this work, we apply a microscopic approach to access the 
time-, momentum-, and angle-resolved dynamics of electrons 
and phonons in optically excited graphene under the influence 
of a variable n- and p-doping. The focus lies in particular on 
the impact of a finite Fermi level on the appearance of the 
technologically relevant carrier multiplication.*^"^® 

This interesting ultrafast phenomenon is related to the lin¬ 
ear electronic band structure of graphene opening up the 
possibility of efficient Coulomb-induced Auger processes. 
A significant multiple carrier generation has been theo¬ 
retically predicted*®’*^’^® and experimentally confirmed in 
graphene.So far, the theoretical studies have been 
constrained to the case of undoped graphene. Introducing a 
non-zero Fermi level in graphene, electrons above the Dirac 
point or above the Fermi level can be considered as charge car¬ 
riers (holes in analogy). In the first case, the carrier multipli¬ 
cation can take place via Auger scattering bridging the valence 
and the conduction band. In the following, we label this pro¬ 
cess as carrier multiplication (CM). On the other side, count¬ 
ing carriers with respect to the Fermi level, the multiplication 
occurs via Coulomb-induced intraband scattering bridging the 
states below and above the Fermi level, cf. Fig. 1. According 
to literature,^® we label this process as hot carrier multipli¬ 
cation (hCM). Here, the actual number of charge carriers re¬ 
mains unchanged in each band. Nevertheless, since these hot 
carriers are crucial for many technological applications, the 
appearance of hCM is also of technological relevance. 


The starting point for the calculation is the many-particle 
Hamilton operator H = Hq-\- H^j + He p + Hc^c, where Hq 
denotes the interaction-free carrier and phonon part, He i the 
carrier-light coupling, Hc.p the carrier-phonon interaction, and 
He c the carrier-carrier interaction.^^ The carrier dynamics is 
described by graphene Bloch equations* corresponding to a 
coupled set of differential equations for the occupation prob¬ 
ability p^{t) in the state k and the band X = (v,c), the mi¬ 
croscopic polarization pk(t) that is a measure for the optical 
transition probability between both bands, and the phonon oc¬ 
cupation nq(t) with the momentum q for different optical and 
acoustic phonon modes 

Pk(f) =(i<+/[ni:''(f)-nr(f)]-ik)pk(f) (1) 

+ '^r(0(Pk(0 ~Pk(0)Fk(t) + ‘^k(t), 

Pk^(f)=±23(nr(f)Pk(f)) (2) 

+r*(',(f)(i-Pk"(f))-r^-;U0Pk(0, 

«q(0 =-7ph(ni(t)-nB) (3) 

+ r-(f)(l+n^(f))-r*^(f)<(f). 

Here, = j(ek~^k) ts the optical transition frequency 
within the linear electronic band structure of graphene 
close to the Dirac point. The carrier-light coupling is de¬ 
termined by (t) = • A(f) with the optical ma¬ 

trix element^^ , the vector potential A(f) represent¬ 
ing the excitation pulse, the free electron mass mo, and the 
charge cq. Here, is the Rabi frequency and 

accounts for intraband transitions.** The many-particle in¬ 
teractions are treated within the second-order Born-Markov 
approximation, which yields a Boltzmann-like scat¬ 

tering equation for the carrier occupation with the time- and 
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FIG. 1. Schematic illustration of Coulomb-induced scattering 
processes for optically excited n-doped graphene. In the presence 
of a finite Fermi level, carrier multiplication (CM) and hot carrier 
multiplication (hCM) need to he distinguished. While the first is in¬ 
duced by Auger processes bridging the valence and the conduction 
band (red arrows) and increasing the number of charge carriers in the 
conduction band, the second corresponds to Coulomb-induced intra¬ 
band scattering (orange arrows) and increases the number of carriers 
above the Fermi level. 


momentum-dependent scattering rates accounting 

for Coulomb- and phonon-induced processes. At the same 
time, the microscopic polarization is damped by the many- 
particle-induced diagonal dephasing Yk{t) and is driven by the 
off-diagonal dephasing term analogy, the equation 

of motion for the phonon occupation nq(f) is obtained and 
contains phonon emission and absorption rates The 

finite phonon lifetime^* 7ph is considered by a coupling to a 
phonon bath at room temperature. The explicit form of the 
time- and momentum-dependent scattering rates is discussed 
in the supplementary material. More details on the diagonal 
and off-diagonal dephasing terms can be found in Malic et 
al.28 


A finite Fermi level Ef breaks the symmetry between the 
valence and conduction band around the Dirac point, cp. 
Fig. 1. As a result, the occupation probability of elec¬ 
trons = p^ and of holes Pk = 1 “ Pk needs to be treated 
separately. As initial condition, we assume a Fermi distri¬ 


bution = 


exp 


{e,h) _ 


where 


k ±EF)/kBT j+l 

stands for the hole and — for the electron occupation at the 
temperature T. Another important aspect of doping is the 
increased screening of the Coulomb interaction. The bare 
Coulomb potential Vq appearing in the Coulomb-induced scat¬ 
tering rates is screened via the dynamic dielectric 

function e(q, co) that is defined by the Lindhard equation,^^’^^ 
cf. the supplementary material for a more detailed discussion. 
Since this many-particle-induced screening is directly influ¬ 
enced by carrier occupations in the conduction and valence 
bands, doping plays a crucial role and has a significant influ¬ 
ence on the ultrafast carrier dynamics in graphene. 


With the presented microscopic approach, we can track 
the relaxation dynamics of non-equilibrium charge carriers in 
time, and energy including the temporal evolution of the car¬ 



FIG. 2. Electron and hole occupations in doped graphene. Angle- 
averaged occupation probability Pk for (a) electrons and (b) holes in 
highly n-doped graphene (Ef = 300 meV) is shown as a function of 
the carrier energy for different times after the optical excitation. Note 
that due to the n-doping the electron relaxation is faster resulting in 
a thermalized hot Fermi distribution already after 30 fs. The inverse 
behavior can be found for p-doped graphene. 


rier density after the optical excitation. First, we focus on the 
electron and hole dynamics in highly doped graphene and then 
we discuss the impact of doping on the appearance of carrier 
multiplication. 


III. ELECTRON AND HOLE DYNAMICS 


Here, we discuss how the doping-induced symmetry break¬ 
ing between electrons in the conduction band and holes in 
the valence band influences the dynamics of optically ex¬ 
cited charge carriers in realistic doped graphene samples. Fig¬ 
ure 2 illustrates the angle-averaged occupation probability 
for (a) electrons and (b) holes for an initial Fermi level of 
300 meV as a function of the carrier energy for different times 
after the optical excitation. Note that for symmetry reasons, 
the physical picture remains the same in p-doped graphene 
samples but electrons and holes switch their roles, respec¬ 
tively. The system is excited by a lOfs pulse with a photon 
energy of 1.5 eV and a pump fluence of 0.3/rJcm^^. The 
characteristics of the excitation pulse correspond to typical 
values that can be realized by standard pulsed lasers.The 
pulse is centered at 0 fs and gives rise to a well pronounced 
non-equilibrium distribution for electrons and holes around 
the carrier energy of 0.75 eV, cf. Fig. 2. For both electrons and 
holes, the efficient carrier-carrier and carrier-phonon scatter¬ 
ing leads to an ultrafast thermalization of the system towards 
a hot Fermi distribution already after some tens of femtosec¬ 
onds. Then, a slower phonon-induced carrier cooling occurs 
that drives the electron and hole occupations towards their ini¬ 
tial thermal Fermi distributions. Due to the increased num- 
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FIG. 3. Doping dependence of carrier multiplication (CM), (a) 

Temporal evolution of the doping-dependent CM for a fixed absorbed 
pump fluence of e^bs = 0.3 /ijcm^^ and an excitation energy of 1.5 
eV. (b) CM as a function of the Fermi level Ef for three fixed time 
delays and (c) the temporal evolution of CM for three fixed Fermi 
levels. Note that CM only takes place for Ef smaller than 0.3 eV and 
for times up to 150 fs. 
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FIG. 4. Auger rates in doped graphene. Temporal evolution of 
Auger rates including impact excitation (IE) and the inverse process 
of Auger recombination (AR) for (a) undoped and (b) highly n-doped 
graphene (e = 300meV). In the latter case, also the corresponding 
intraband Auger rates (Jie, )^«) describing Coulomb-induced scat¬ 
tering bridging the states below and above the Fermi level are illus¬ 
trated (dashed lines), cf. Fig. 1. In undoped graphene, there is a 
clear asymmetry between both Auger processes in favor of IE in the 
first 100fs. In the doped case, is significantly larger than //£, 
however, for intraband Auger rates, the situation is opposite and //£ 
clearly prevails over resulting in a pronounced hot carrier multi¬ 
plication. 


ber of available scattering partners in the conduction band of 
n-doped graphene, the Coulomb-driven carrier thermalization 
occurs faster for electrons. Here, a hot thermalized Fermi dis¬ 
tribution is already reached after 30 fs, while at the same time 
the holes exhibit still a non-equilibrium distribution, cf. the 
purple lines in Fig. 2. 


IV. CARRIER MULTIPLICATION (CM) 

Now, we study the impact of an initial Fermi level fs/r on the 
Coulomb-induced multiple carrier generation, which is gener¬ 
ally defined as the ratio between the number of overall gener¬ 
ated electron-hole pairs and the optically excited charge carri¬ 
ers 

{h)CM='^^, (4) 

^opt 

where n is the total carrier density, tif the initial thermal car¬ 
rier background, and tiopt the optically excited carrier den¬ 
sity. All contributions contain both electrons in the conduc¬ 
tion band as well as holes in the valence band. For doped 
graphene, the specific definition depends on the physical situ¬ 
ation: For optical measurements probing vertical carrier tran¬ 
sitions, a definition with respect to to the Dirac point is rea¬ 
sonable, i.e. = p^ and pj^' = 1 — p^’. The carrier density 


then reads 

E Pt ( 5 ) 

X=£,,h‘, k 

where is the graphene area and CL (dy) denotes the spin 
(valley) degeneracy. On the other side, electric transport phe¬ 
nomena are also of great interest, where hot electrons around 
the Fermi level are relevant giving rise to a hot carrier multi¬ 
plication, cf. Fig. 1. This situation will be discussed in the 
next section. 

Treating the full set of graphene Bloch equations, we have 
microscopic access to the temporal evolution of the carrier 
density including the contributions of carrier-light, carrier- 
carrier, and carrier-phonon interactions. Figure 3(a) illustrates 
the temporal evolution of CM as a function of the initial Fermi 
level Ef at a fixed absorbed pump fluence of 0.3 pJcm^^. The 
surface plot reveals that doping clearly reduces the CM effi¬ 
ciency: The lower Ef, the higher is the CM factor reaching 
values of up to approximately 1.7 for undoped graphene at the 
considered pump fluence, cf. Fig. 3(b). CM can be observed 
for Ef of up to 300 meV. Its maximal lifetime is approx¬ 
imately 150 fs for undoped graphene and becomes signifi¬ 
cantly shorter for increasing doping, as illustrated in Fig. 3(c). 
The observed CM in the low-doping case can be explained 
by the strongly efficient impact excitation (IE) prevailing over 
the inverse process of Auger recombination (AR), which is a 
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FIG. 5. Doping dependence of hot carrier multiplication (hCM). 

The same plot as in Fig. 4 illustrating now the temporal evolu¬ 
tion of the doping-dependent hCM for a fixed absorbed pump flu- 
ence of e^bs = 0.3 /tJcm^^. The inset in (c) shows a direct com¬ 
parison between theoretically predicted (solid lines) and experimen¬ 
tally measured time-resolved ARPES data (dashed lines) for: (i) n- 
doped graphene {Ef = 380 meV, Cabs =0.5 /t Jcm^^) and (ii) p-doped 
graphene {Ep = 240 meV, Eabs = 1.5/iJcm^^). The experimental 
data is taken from Johannsen et al.*^. 


result of the large gradient in carrier occupation around the 
Dirac point, cf. Fig. 2(a). For undoped graphene, the prob¬ 
ability for IE can be written as p’’(l — p'^) « 1, whereas the 
probability for AR is given by p'^(l — p'’) « 0 (since p'' « 1 
and p'^ « 0). In this case, IE is significantly favored by Pauli 
blocking during the initial dynamics. This is reflected by the 
corresponding rates Jie and that are shown in Eig. 4(a). 
We observe that the IE rate is clearly higher for a time range 
of approximately lOOfs determining the strength of the ap¬ 
pearing CM. During the carrier relaxation both rates converge 
to the same value and end up in an equilibrium, where no 
more carriers are generated. The lifetime of the CM is deter¬ 
mined by the duration of the imbalance between IE and AR 
rates in combination with the interplay with competing chan¬ 
nels of carrier-phonon scattering, which transfer energy from 
the electronic system to the lattice. With an increasing dop¬ 
ing, Auger scattering becomes more and more Pauli blocked 
resulting in overall lower rates, cf Pig. 4(b). PorPermi levels 
higher than 300 meV, AR becomes the dominant relaxation 
channels and CM does not appear anymore. 


V. HOT CARRIER MULTIPLICATION (HCM) 

Now, we focus on the situation, where charge carriers are 
defined with respect to the Permi level, i.e. for n-doped 
graphene the upper Dirac cone is split into p^’^ = p^ for k> kp 
and p''’'^ = l-p^fork < kp with the Permi momentum kp. 
The bottom cone remains unaffected with p^ '’ = p^ — 1— p^. 
The according hot carrier density h is given by 

_ _ / k<kp k>kf \ 

(6) 

and the associated hot carrier multiplication is highly relevant 
for transport phenomena, cf. Pig. 1. Note that we ob¬ 
tain symmetric results for n- and p-doped graphene, since the 
contribution of both electrons and holes to the carrier density 
is considered. Por undoped graphene, both definitions of car¬ 
rier density [Eqs. (5) and (6)] and carrier multiplication are 
equivalent for symmetry reasons. 

The surface plot in Pig. 5(a) illustrates hCM as a func¬ 
tion of the Permi level Ep and time at a fixed absorbed pump 
fluence of Cabs = O.SpJcm^^. In contrast to the behavior of 
CM, we observe a clear increase of hCM with doping. We 
reach hCM factors of up to approximately 2 (at the consid¬ 
ered pump fluence) with a lifetime of about 200 fs for highly 
doped graphene with Ep = 300 meV, cf. Pigs. 5(b) and (c). 
There is nearly a linear dependence between hCM and doping: 
The smaller Ep, the less pronounced is hCM, and the shorter 
is its lifetime. The probability for intraband IE processes is 
given by p^'’‘^(l — which is initially large compared to 
the probability for intraband AR processes p'^’‘^(l — p*’'^), cf. 
Pig. 1. With the increasing Permi level, the intraband Auger 
processes are shifted into a region of higher density of states 
making them more efficient, as reflected by the much higher 
scattering rates displayed in Pig. 4(b). The initial strong im¬ 
balance between Yar and fiE gives rise to a pronounced hCM. 

Besides the discussed doping dependence, CM or hCM are 
strongly sensitive to the excitation regime. A detailed dis¬ 
cussion is provided in the supplementary material, where a 
semi-analytical approach is presented focusing on the purely 
Coulomb-induced CM and hCM. 


VI. DIRECT COMPARISON TO EXPERIMENTAL DATA 

After having presented the theoretical results on the doping 
dependence of the carrier multiplication, we perform a direct 
comparison with recently performed time-resolved ARPES 
measurements on n- and p-doped graphene samples.*^ We 
explicitly take into account the experimental conditions, such 
as the Permi level and the excitation strength. Considering 
that our microscopic theory does not contain any fitting 
parameters, we obtain an excellent agreement between theory 
and experiment, cf. the inset in Pig. 5(c). There is a clearly 
higher hCM for n-doped graphene reaching values of up 
to 2.2 in the theory and more than 3 in the experiment. In 
contrast, for p-doped graphene only a small hCM of 1.4 or 
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1.2 is obtained in theory and experiment, respectively. This 
pronounced difference is not due to the type of doping (n, 
p), as one might assume.*^ It can be clearly explained by 
the differences in the applied fluence Cabs and the actual 
Fermi level Ep. Note that at the exactly same conditions 
with respect to Ep and Cabs, we obtain the same hot carrier 
multiplication for both n- and p-doped samples. However, 
the experiment has been performed for: (i) n-doped graphene 
with the Fermi level Ep = 380 meV and an absorbed pump 
fluence of eab,s = 0.5 /r Jcm^^ and (ii) p-doped graphene with 
Ep = 240 meV and gabs = 1.5/rJcm^^. As shown in Fig. 
5(b), hot carrier multiplication increases almost linearly with 
the Fermi level. Furthermore, it is strongly suppressed in the 
strong excitation regime, i.e. the larger the pump fluence, 
the less efficient is the hCM, as illustrated in Fig. S3 in the 
supplementary material. As a result, the n-doped graphene 
sample shows a much more pronounced hCM, since its Fermi 
level Ep is significantly higher and since the experiment has 
been performed at a clearly smaller pump fluence compared 
to the p-doped graphene sample. 

In summary, we have presented a microscopic study of the 
carrier dynamics in doped graphene samples, in particular fo¬ 
cusing on the impact of a finite Fermi level on the (hot) carrier 
multiplication. We reveal the appearance of Auger-induced 
carrier multiplication up to Fermi levels of 300 meV. In the 
case of the hot carrier multiplication occurring within one 
band doping is even advantageous, since it increases the phase 
space by providing a large number of available scattering part¬ 
ners. Finally, we have directly compared our results to re¬ 
cent time-resolved ARPES measurements finding an excellent 
agreement and providing a microscopic explanation for the 
observed different behavior in n- and p-doped graphene sam¬ 
ples. Our results contribute to a better understanding of the 
ultrafast carrier dynamics in realistic graphene samples and 
give valuable insights into the technologically relevant carrier 
multiplication in graphene. 
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FIG. 6. Analytically obtained CM and hCM as a function of doping 
at a fixed pump fluence (solid line) of e^bs = 0.1 p Jcm^^ compared 
to the numerical solution (dots) obtained by evaluating the graphene 
Bloch equations. 


analogous way. Our approach is to evaluate the total carrier 
density n and the energy density S' of the electronic system: 


L Pk and <r= 


A=e,/i;k 


L2 


(7) 


/tk 


Assuming a Fermi distribution with Ep ^0 for the carrier oc¬ 
cupation p^, we obtain the carrier and energy density as a 
function of the temperature T and the Fermi level Ep: 

OiiOx) w—1 k 


2k 


dk- 


(7=± ' 


0 g{kvF+aEF)l(kBT) _|_ l 


TTV 


o'=± 


( 8 ) 


SiT,Ep) = 


GsOvVp 

2 % 

Ak 




dk 


g{kvF+aeF)l{kBT) _|_ \ 




o'=± 


(9) 
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APPENDIX 

1. Semi-analytic approach to purely Coulomb-induced CM 

To be able to find optimal excitation conditions and the dop¬ 
ing regime for (hot) carrier multiplication, we present a semi- 
analytical approach considering the purely Coulomb-induced 
carrier dynamics within statistical methods. 

Here, we focus on the carrier multiplication (CM), but the 
processes of hot carrier multiplication can be treated in an 


where Li„{.) denotes the polylogarithm function. 

After an optical excitation, the electronic system is de¬ 
scribed by the carrier density Uexc = no -t- nopt and the energy 
density Sgxc = Sq + ^abs- Here, ^abs denotes the absorbed 
pump fluence and riopt = 2(?ab.s/^® the optically excited car¬ 
rier density with the photon energy hco. Eurthermore, the 
thermal densities no = n{To,Ep) and Sq = S{Tq,Ep) are de¬ 
termined by the initial temperature To and Eermi level Ep. To 
access the CM factor we need to know the final state of the 
carrier system after the Coulomb-induced thermalization. We 
exploit the fact that the Coulomb dynamics conserves the to¬ 
tal energy density. Due to Auger processes, the total carrier 
density is not conserved, but instead the Eermi level. Using 
Eq. (9), we are able extract the final temperature 7y of the hot 
thermalized system by numerical variation with the condition 
that S{Tf,Ep) = Sexc- Having the final temperature, we can 
also access the final carrier density rif = n{Tf,Ep) and there¬ 
with also the CM factor, cf. Eq. (4). Eor undoped graphene. 
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FIG. 7. (a) Carrier multiplication (CM) as a function of Fermi level 
and pump fluence at room temperature and at an excitation energy 
of 1.5 eV. (b) Doping-dependent CM for three representative pump 
fluences and (c) fluence-dependent CM for three fixed Fermi levels. 


the purely Coulomb-induced Tf, nf, and CM can be obtained 
analytically.'^ Note however that the approach does not in¬ 
clude carrier-phonon scattering, which cannot be treated stat¬ 
ically. Nevertheless, Fig. 6 shows an excellent agreement 
between the analytically [cf. Eqs. (8)-(9)] and numerically 
[cf. Eqs. (l)-(3)] obtained doping-dependent carrier multi¬ 
plication at a fixed pump fluence of Cabs =0.1 /f Jcm^^. The 
higher the Eermi level, the smaller is CM, and the larger is 
hCM. Eor highly doped graphene samples with £p > 0.4eV, 
CM becomes smaller than 1. In contrast, hCM reaches val¬ 
ues of almost 5. The same qualitative dependence on dop¬ 
ing is found by numerically evaluating the graphene Bloch 
equations (shown in the main part), which also include the 
carrier-phonon coupling. This demonstrates that the presented 
semi-analytical approach already sufficiently covers the most 
important aspects of doping-dependent carrier multiplication. 
Eor quantitative insights and time-dependent CM, the full mi¬ 
croscopic approach is required. 


2. Fluence dependence of carrier multiplication 

Now, based on the semi-analytical approach presented in 
the last section, we discuss the impact of the excitation 
strength on the appearance of the CM and hCM. Eigures 
7(a) and 8(a) show surface plots illustrating doping and flu¬ 
ence dependence of CM and hCM at 300 K and at an excita¬ 
tion energy of 1.5 eV, respectively. Our calculations demon¬ 
strate that both CM and hCM are strongly suppressed in the 
strong excitation regime. Eor pump fluences higher than ap- 


FIG. 8. The same as in Fig. 7, however illustrating doping and flu¬ 
ence dependence of the hot carrier multiplication (hCM). 

proximately 5/iJcm^^ (hot) carrier multiplication does not 
occur any more. The larger the pump fluence, the more 
scattering partners are available giving rise to a faster car¬ 
rier dynamics.'® As a result, the thermalized distribution is 
reached within the first tens of fs and the asymmetry of im¬ 
pact excitation and Auger recombination rates vanishes much 
faster. As a result, carrier multiplication clearly decreases at 
enhanced pump fluences, cf. Eigs. 7(c) and 8(c) Eurthermore, 
Eigs. 7(b) and 8(b) also illustrate the decrease (increase) of 
CM (hCM) with the Eermi level, as already discussed in the 
main part. 


3. Dynamical screening 


The efficiency of Auger processes is sensitive to the dynam¬ 
ical screening of the Coulomb interaction that is expressed by 
the Lindhard equation:^^ 


e(q,a)) = l-2V^Y^ a' 

AA'k ^k+q - Ek + + '4,q 


-Pk 


( 10 ) 


where O) corresponds to the energy transfer of the scattering 
processes and the factor E^^^ = |/£/^r'F^*(r)e'‘' ’’'F^lq(r)|^ 
includes the tight-binding wave functions (r). 

The appearing momentum-dependent dephasing can be 
derived from the dynamics of inhomogeneous charge fluctu¬ 
ations (Pq) = -^Lk(ak-q“k) = -^LkC7k-q,k,- The Corre¬ 
sponding equation of motion reads:^^ 


J<7k-q,k — * (£k-q £k)<7k-q,k + * (Pk Pk-q) *^k'-q,k'■ 
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Using the ansatz 0'k-q,k(f) = e^'(®+'^^'ak-q,k(0), a solu¬ 
tion for <7k-q,k can be obtained and leads to the dynam¬ 
ical screening expressed by the Lindhard formula with a 
fluence-dependent dephasing To obtain the momentum- 

dependent we write the in-coherent equation of motion 

for the electron charge density in second-order Born-Markov 
approximation; 

d 7t 

;^0-k-q,kInc =-T E i^D + £c - £b - £a) X (H) 
at n 

bed 

abd 

This equation contains the momentum-dependent dephasing 
term, which is given by ^(Jk-q,k|ne 4.qO-k-q,k- Neglecting 
all polarization terms Oab —> Pa, the dephasing 5k.q yields: 

4^q = 5”;f+5™f (12) 


with the in- and out-scattering rates 


q E^k,q,a,b,c (1 ~ Pa)PbPc5(ek TEa + £b + £c)) 

" abc 

‘^k'q’"'^ ~T E^k>q’a.b,cPa(l — Pb)(l ~ Pc)8{e^ -f Ca -f £b + £c) 

" abc 


with Vk n a b C = t that 

■4? ? . k A k('A(' k^A^k^A^ k qAkcAc 

also directly appear in the graphene Bloch equations, cf. Eq. 
(2). Calculating these scattering rates, we have also a consis¬ 
tent access to the dephasing term influencing the dynamical 
screening of the Coulomb interaction. Note that a vanishing 
dephasing would suppress the parallel Auger scattering con¬ 
tributions in graphene due to a diverging e(q, O)).^* However, 
in realistic cases the many-particle processes lead to a signifi¬ 
cant dephasing opening up the Auger channels for the carrier 
relaxation. 
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